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Abstract 



We summarize our recent results on the large N renormalization group 
(RG) approach to matrix models for discretized two-dimensional quan- 
tum gravity. We derive exact RG equations by solving the reparametriza- 
tion identities, which reduce infinitely many induced interactions to a fi- 
nite number of them. We find a nonlinear RG equation and an algorithm 
to obtain the fixed points and the scaling exponents. They reproduce the 
spectrum of relevant operators in the exact solution. The RG flow is vi- 
sualized by the linear approximation. 
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Abstract 

We summarize our recent results on the large N renormalization group (RG) 
approach to matrix models for discretized two-dimensional quantum gravity. 
We derive exact RG equations by solving the reparametrization identities, 
which reduce infinitely many induced interactions to a finite number of them. 
We find a nonlinear RG equation and an algorithm to obtain the fixed points 
and the scaling exponents. They reproduce the spectrum of relevant operators 
in the exact solution. The RG flow is visualized by the linear approximation. 

1. Introduction 

Matrix models have been a powerful tool to study two-dimensional quantum gravity (2D QG) 
via simplicial decomposition of the spacetime 1 . We are particularly interested in the 2D QG 
both as a string theory and a toy model for the QG in higher dimensions. Exact solutions of the 
matrix model have been successfully obtained for 2D QG coupled to minimal conformal matter 
with central charge c < 1, but unsuccessful for c > 1 cases which correspond to realistic string 
theories. We need to obtain approximation schemes which provide us with correct results for the 
exactly solved cases and enable us to calculate critical coupling constants and critical exponents 
for unsolved matrix models, especially for c > 1. 

A large A renormalization group 2 (RG) has been proposed by Brezin and Zinn- Justin as such 
an approximation method 3 , identifying the size of the matrix with the inverse lattice constant 
of the simplices, and was discussed by several groups 4 . We have succeeded to derive explicitly 
exact RG equations for 0(A)-vector models 5 and one- 6 and two- matrix models 7 . We found that 
it is crucial to take account of the reparametrization identities in order to obtain a meaningful 
RG equation. The RG equation for the matrix model turned out to be nonlinear in contrast to 
the linear RG equation for the vector model. Moreover, we found that the global picture of the 
RG flow can be drawn practically by a linear approximation to the nonlinear RG equation 7 . 

2. Nonlinear RG equation for the one-matrix model 

The partition function Z N (gj) of a one-matrix model with a general potential V(<p) = J2k>i 
is defined by an integral over an A x A hermitian matrix 0, Z^(gj) = f d<p e~ NtrV( -^\ Starting 
from an (iV+1) x (A + l) hermitian matrix variable $, we decompose it into an A x A hermitian 
matrix 0, a complex A-vector v and a real scalar a. Taking cubic interaction V((f>) = + |0 3 , 
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we integrate over v exactly to obtain 

Z N+1 ( g) = (_L_y J d(f) e -(N+l)trV(4>) . J d(y e -(iV + l)V(a)-trlog(l+^ + a))_ (1) 

We can evaluate the a-integral by the saddle point method, systematically in 1/iV-expansion. If 
we explicitly retain only the leading part of the 1/iV-expansion, we obtain a difference equation 
relating Z N+1 and Z N . We shall denote the average Z^ 1 J d<p{- ■ ^e~ NtrV ^ by (■ ■ •). Factorization 
holds in the large- N limit (O O') = (0)(0') + 0(N- 2 ) for a multi-point function of ^(AO- 
invariants O, O'. By using this property we obtain for the free energy F(N,g) = 
-N~ 2 log(Z N (g)/Z N (0)) 



F(N, g ) = -l + (±tr V(<j>)\ + V((a s )) + / i tr log (1 + g(cf> + (a a }))\ + O (i 



(2) 

Here a s = a s (<p) is determined as a [/(iV)-invariant function by the saddle point equation 

N (l/g + a B ) +<j) 

The right hand side of eq. (fj) consists of products of (-^ trft) = jdF/dgj, j — 1, 2, . . .. The 
crucial observation we make is that there is an ambiguity to define a flow in the coupling constant 
space, since we can make an arbitrary reparametrization of matrix variable <ft in the partition 
function. This will enable us to express the <9-F/<9g,,-dependence (j ^ 3) in terms of dF/dg^, and 
to reduce all the induced interactions to those in the original potential. The resulting identities 
are a part of the discrete Schwinger-Dyson equation for the system. 

More concretely, if we perform reparametrizations regular at the origin, <$' = <fi + e<f) n+1 (n > 
—1) for the partition function, we obtain identities 

? \iv tr0J iv tr<r ~ 7 = \iv tr (> n+v ^))) (» ^ ( 4 ) 

These identities imply relationship among dF/dgj's in the large N limit, 

« <^ ^ •/ . « dF OF <9F 

2rl ^— + Z^^( n_ ^)^— a = 2J n + %^ ■ ( 5 ) 

og n J^l dgjdgn-j ~ d5f n+fc 

By introducing the resolvent W^fV) = (l/iV)tr(,2 — </>) _1 , eq. @ is neatly summarized into a 
single equality called the loop equation 1 

W(z)) 2 - V'(z) (w(z)) + Q(z) = 0, Q(z) = E Vik+ ^ Z) (± tr (0 - z)^ 1 ) . (6) 

Eqs. hold for a generic polynomial potential V(<p) = J2T=i ^<^ fc > f° r our case m = 3, 

is given by Q(z) = 1 + gz — g 2 + 3g 3 dF/dg. Consequently we obtain a differential equation 
obeyed by the free energy F(N, g) of the one-matrix model with the cubic coupling, 



n Wm + 2 



F(N, g) = G (g, ^\ + O (J^j , a(g, a) = (a s ) = -g + 3g\ 

—l/g-a 



G(g,a) = -?-a + -a 2 + ^a 3 + log (I + ga) + / dz ({W(z, g, a)) - - . (7) 



To establish an algorithm for fixed points and scaling exponents, we first concentrate on the 
leading part F°(g) of the free energy in the 1/iV-expansion. It is easy to see that F° satisfies 
2F°(g) = G (g,dF°(g)/dg). We assume that F° consists of regular and singular parts around 
a fixed point F°(g) = E^o^G? - 9*) k + Ek=oh(g - g*) k+2 ^° (2 - 7 o t N), and that 
G(g,a) is regular around (g*,ai). 7 o is referred to as the susceptibility exponent. By comparing 
the coefficients of various powers of g — g* on both sides of the RG equation, we have a set 
of equations to determine unknowns 70, g*, and bk- We can also derive the RG equation 
for higher genus contributions and find that they exhibit the double scaling behavior under a 
plausible assumption. 

We can as well use the eigenvalue representation for the RG equation which is particularly 
useful for multi-coupling cases. We find four fixed points (Yang-Lee edge singularity, pure gravity 
and Gaussian) in the two-coupling (g 3 and g 4 ) case 6 ; there the spectrum of relevant operators of 
the exact solution (70 = —1/3,-3/2 for Yang-Lee and 7o = —1/2 for pure gravity) are exactly 
reproduced. 



3. The two-matrix model 



+ |0 3 is defined as 



The partition function for the two-matrix model with V((j),g) - ... 

Z N (g+,g-,c) = y# + d0_exp[-7Vtr(y(0 + , 5+ )+y(0_,^_)+c0 + 0_)]. (8) 

Similarly to the one-matrix case, we integrate the N + 1-th row and column exactly and evaluate 
the singlet integral by the saddle point method. We consider the case g + = g- = g which 
describes the Ising model coupled to 2D QG in the absence of an external magnetic field. We 
assume that (tr0+) = (tr0i) holds for any j, and solve the reparametrization identities (a 

fourth-order algebraic equation for the resolvent in need, Wo(z) = ^Tr^ Z ~^^ + ^ A 

Then we find a nonlinear RG equation for the two-matrix model 



z + 



d 

N w +2 



G(g, c; a g , a c ) 



/ dF dF\ /l 
-^a g + (1 + c)a 2 + ^-a 3 + 2 log(l + go) 



9 

tun 

2 9 

rl/g+a / / - , .\ 2 

+ / dz (W (z;g,c;a g ,a c ) 



a(g,c;a g ,a c ) = (a± >s ) 



9 



1 + c 



-1 + c a c + 



3g 



(9) 



We find three fixed points (critical Ising, pure gravity and Gaussian) in the physical region of 
the coupling constant space 7 ; there the spectrum of relevant operators of the exact solution 
( 7 o = —1/3, —5 for critical Ising and 7 o = —1/2 for pure gravity) are again exactly reproduced. 



4. Linearized RG flow 

We observe generically that the coefficients of the nonlinear terms in the G-function, d n G/da n (g, 0) 
for n > 2, are suppressed by increasing powers of g. Therefore we expect that within a region 
g < 1 where most fixed points lie, the RG flow should be well approximated by the linearized 
/5-function, (3n n (g) = dG/da(g,0). Here we exhibit the RG flow in one- (Fig.l) and two-matrix 
models (Fig. 2) described by the real parts of (3\i n 7 . We find that the critical lines (real lines) ema- 
nating from the multicritical fixed Doints toward the Dure gravity fixed Doint are characterized as 



the RG trajectories corresponding to the least relevant perturbations of the UV theories. These 
flows are consistent with those for matter fields over a fixed background, and suggest a gravita- 
tional analogue of the celebrated c-theorem 8 stating the decrement of the degrees of freedom of 
the matter sector along the RG flow. 
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Figure 1: RG flow for the one-matrix model. Figure 2: RG flow for the two-matrix model. 

The arrows represent a flow from the UV (N = oo) to the IR (N = 0). 



